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Abstract 

The scale invariance of the 0(3) sigma model can be broken by gauging a U{1) subgroup of the 
0(3) symmetry and including a Maxwell term for the gauge field in the Lagrangian. Adding also a 
suitable potential one obtains a field theory of Bogomol'nyi type with topological solitons. These 
solitons are stable against rescaling and carry magnetic flux which can take arbitrary values in 
some finite interval. The soliton mass is independent of the flux, but the soliton size depends on it. 
However, dynamically changing the flux requires infinite energy, so the flux, and hence the soliton 
size, remains constant during time evolution. 
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1 Introduction 



The 0(3) sigma model in (2+1) dimensions is a popular model in theoretical physics. Stat- 
ically it is integrable and of Bogomol'nyi type, i.e. all minimal energy solutions can be 
obtained by solving the first order Bogomol'nyi equations (which imply the second order 
Euler-Lagrange equations). As a result one can explicitly write down soliton solutions of 
arbitrary degree in term of rational functions From the point of view of a particle physi- 
cist, however, the model has one important drawback: it is scale invariant and as a result 
its soliton solutions have an arbitrary size, making them unsuitable as models for particles. 
Numerical simulations of the solitons' interaction behaviour in the (2+l)-dimensional model 
suggest that the solitons do indeed change their size during interactions and generically turn 
into singular configurations of zero size p[ . The most obvious way to break the scale invari- 
ance of the model is to add terms to the Lagrangian density which contain a different number 
of derivatives from the sigma model term (which contains two). Indeed, the inclusion of a 
Skyrme term (four derivatives) and a potential term (no derivatives) leads to so-called baby 
Skyrme models which have soliton solutions of definite size. Such models are close analogues 
of the (3+1) dimensional Skyrme model and therefore physically interesting, but they are 
neither integrable nor of Bogomol'nyi type and can only be studied with considerable nu- 
merical effort IP] . A mathematically more elegant way of breaking the scale invariance of the 
0(3) sigma model is to add a potential term only and prevent the solitons from collapsing 
by making them spin. In [^] it was shown that with a suitable choice of the potential such 
a model is of Bogomol'nyi type. Its soliton solutions, called Q-lumps, can be written down 
explicitly and their interaction behaviour was studied in 0]. 

Here we investigate the possibility of breaking the scale invariance of the sigma model 
by introducing a U{1) gauge field whose dynamics is governed by a Maxwell term. This 
possibility is also potentially of interest in the (3+l)-dimensional Skyrme model. There the 
gauging of a U{1) subgroup and the inclusion of a Maxwell term is physically natural, and it 
would be aesthetically appealing if one could do away with the Skyrme term and still retain 
a model with stable soliton solutions. Since in three spatial dimensions the sigma model 
term has the opposite scaling behaviour from the Maxwell term, simple scaling arguments 
do not rule out solitons of a definite size in a "Skyrme-Maxwell model without a Skyrme 
term". At the end of the paper we will make some conjectures about this possibility. 

Since the solitons here carry magnetic flux it is appropriate to compare them to vortices. 
Vortices either have quantised flux in which case they are topologically stable (e.g. in the 
abelian Higgs model, see and references therein) or they have arbitrary flux in which case 
they are not topologically stable (e.g. non-topological Chern-Simons vortices, see 0). The 
topological stability of the solitons studied here, however, is independent of their magnetic 
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flux. Thus they can have arbitrary flux and yet be topologically stable. 



2 Topological stability 

In the 0(3) sigma model the basic field is a map from (2+l)-dimensional Minkowski space 
to the 2-sphere of unit radius. Here Minkwoski space is assumed to have the signature 
(—,+,+), and its elements are written as {t,x) or alternatively x", a = 0,1,2; for partial 
derivatives with respect to these coordinates we write da- The field, denoted <p, has three 
components 0i, 02 and 03 satisfying the constraint <p ■ 4> = (pf + (p"^ + (pl = 1- The potential 
energy functional is 

EM = \j d'x {i^^cj,f + id2cf>f) . (2.1) 
To obtain finite-energy configurations one requires that, at all times t, 

lim cb(t,x) = n, (2.2) 

where n is a constant unit vector which we take to be n = (0, 0, 1) for definiteness. This 
condition allows one to add the point oo to physical space R^, thus compactifying it to a 
topological 2-sphere. As a result a field at a fixed time may be viewed as a map from 
one 2-sphere to another and therefore has an associated degree deg[0]. This degree is a 
homotopy invariant and therefore cannot change during time evolution. 

For our purposes it is convenient to express the degree in terms of the current 

A;a = ^ea/3-,0-5>x9^0, (2.3) 

where e^/?^ is the totally antisymmetric tensor in three indices, normalised so that eoi2 = 1- 
One finds 

deg[0] = ^ fd^xko = ^ [ (fx (j)-di(j) X ^20. (2.4) 
It is easy to see that the divergence of ka vanishes independently of the equations of motion; 



together with the boundary condition ( p.2| ) this explicitly shows the conservation of the 
degree. The degree of a configuration is also important because it provides a lower bound 
on its energy |^: 

^.[0] >47r|deg[(/)]|. (2.5) 



The energy functional E„ and the boundary condition ( p.2| ) are invariant under the group 



of global rotations of the field cf) about the fixed vector n. This is the U{1) symmetry we 
want to gauge. Thus we introduce a f/(l) gauge field and a covariant derivative 

Da(t> = da(f> + AaTl X (f). (2.6) 
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Defining the field strength as usual via Fais = daAp — dpA^ with magnetic field F12, we can 
write down the potential energy functional which is the subject of this paper: 

= \jd\ {{D,ct>f + (^20)' + (1 - n-(/>)2 + F^J . (2.7) 

In the gauged model the topological current ka^ while still divergence free, is unsatisfac- 
tory because it is not gauge invariant. Through trial and error one finds that the current 

= \e.p, {(t>-D^(t> X D'<ct> + ^'^^(l - n-0)) , (2.8) 

which is manifestly gauge invariant, differs from only by the curl of another vector field 

3a = k^ + e^p^d^{{l-u-4>)A~<). (2.9) 

Thus ja also has vanishing divergence, daj°' = 0, and together with the boundary condition 
(|2.2[) this implies that the degree of a configuration cf) can be expressed as 



deg[0] = ^ / d'^x]^ = 1-Jd^x {(t>-Di4> X D2ct) + ^12(1 - n-(^)) . (2.10) 
Returning to the energy functional -Egaugo we note that it can be re-written as 

Ega,gJ(/>,A,A2] = \j ((Di0±0x D20)' + (Fi2T(l-n-<^))2) 

± j £x {4>-D^cj> X D^cf) + Fi2(l - n-0)) , (2.11) 
where we used that (f)-Da(f> = 0. Together with the formula for the degree ( |2.1CI| ) this implies 

Ega,g40,Ai,A2] >47r|deg[(^]|, (2.12) 
with equality if and only if one of the Bogomol'nyi equations holds: 

Di(f) = IF0 X L)20 

F12 = ±(l-n-0). (2.13) 

It is instructive to write these equations also in a different form which results when the 
target space S"^ is stereographically projected onto C U{oo}. More precisely defining a 
complex-valued field u = Ui + iu2 via 

1 + 03 1 + 03 

the Eqs. ( |2.13| ) become, in terms of u, 

Diu = ^iD2U 

F12 = ±7^, (2.15) 
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where Dj now stands for dj + iAj, j = 1, 2. In the gauge diAi + ^2^2 = these equations 
imply the following second order equation for = Inu: 

Ay? = ^r^- (2-16) 

The integrability of such "non-linear Laplace equations" has been studied in the literature, 
but the present equation lies outside a small class of such equations which are known to be 
integrable by standard methods, such as scattering transforms @]. To find solutions of the 
Bogomol'nyi equations ( p.l3| ) we therefore resort to numerical methods. 



3 Solving the Bogomol'nyi equations 



When seeking solutions of Eq. ( p. 131) with non-zero degree we restrict attention to fields 



which are invariant under simultaneous rotations and reflections in space and target space. 
Thus we assume that (p is of the so-called hedgehog form 

(l){x) = {sin f{r) cos NO, sin f{r) sin N6, cos f{r)), (3.1) 

where (r, 6) are polar coordinates in the a;-plane, is a non-zero integer and / is a function 
satisfying certain boundary conditions to be specified below. The gauge field is assumed to 
have only a 6'-component which is of the form 

Ae = Na{r). (3.2) 

(For a more detailed justification of this ansatz see ||^.) Then the magnetic component of 
the field strength is simply 

F12 = N-. (3.3) 

r 

To obtain fields which are regular at the origin we require 

/(O) = 71 and a(0) = 0, (3.4) 
and to ensure also that the energy is finite we impose 

lim f(r) = and lim a(r) = 0. (3.5) 

One checks that the degree of such a configuration is —N. 

The Bogomol'nyi equations ( p. 131 ) imply the following coupled first order differential 
equation for / and a: 

f = _|iv|^sin/ (3.6) 
r 

T 
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where the alternative signs of ( p.l3| ) have been absorbed into the modulus sign. For brevity 
we will refer to the boundary value problem posed by these differential equations together 
with the boundary conditions ( p.4| ) and ( |3.5| ) as BVP. As a first step in its discussion we 
establish the 

Proposition: The boundary value problem BVP has no solution if \N\ = 1, but it has a 

one-parameter family of solutions if\N\ > 1. 

Before entering the proof we note the solutions of Eqs. (|3.6| ) and (|3.7| ) for small r. Using 
the boundary condition (|3.4|) and keeping only the leading powers in r one finds 



/^TT + Arl^l and a ^ -^r^, (3.8) 

where A is an arbitrary constant. When integrating Eqs. ( |3.6| ) and ( p.7| ) numerically one 
cannot start the integration at the regular-singular point r = 0. Instead we integrate from 
some small value (r = 10"^ in practice) outwards, imposing the initial values there according 
to ( |3.8| ). According to the proposition, there is a family of values for A that will lead to a 
solution satisfying the boundary conditions p.5| ) at infinity when |A^| > 1 but there is no 
such value when |A^| = 1. 

The proof of the proposition proceeds in four steps; in it a function g is called increasing 
(decreasing) ii x > y ^ g{x) > g{y) {x > y ^ g{x) < g{y))- 

1. It is clear from Eq. (|3.7|) that a' < 0, i.e. a is a decreasing function. Moreover, since 
< (1 — cos/) < 2, the boundary conditions (|3.4|) and the intermediate value theorem imply 
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2|Ar| 



< a < 0. (3.9) 



2. For any solution of BVP, /(r) G (0, vr] for all r. To prove this we first show that /(r) > 
for all r. For suppose that / were close to zero, so that Eqs. (|3.6| ) and (|3^ ) simphfy to 



/' = -|JV|^/ (3.10) 

a' = -^f-. (3.11) 

Then, from the first equation f d\nf = —\N\ J dr {a + l)/r. However, while the left hand 
side diverges as / 0, the right hand side is finite for any finite interval of integration (we 
can assume without loss of generality that r = is not included in the interval). Thus / 
cannot vanish for any finite value of r; since it is continuous and non-zero for r = it is 
positive for all r. By a similar argument one shows that if /(^q) > vr for some rg > then 
/(r) > vr for all r > tq, which violates the boundary condition ( p.5|) . Thus < / < tt as 
claimed. 

3. For any solution of BVP, a(r) > —1 for all r. Suppose this were not the case. Then 
there is an ri such that a(ri) < —1 and, since a is decreasing, a(r) < —1 for all r > ri. 
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However, since / G (0, tt] , Eq. ( p.6| ) then implies that / is increasing for r > ri, which is 
incompatible with the boundary condition ( |3.5| ). Combining this result with the inequality 
( p.9| ) we conclude — 1 < a < 0. Note that it then follows from Eq. ( |3.6| ) that / is a decreasing 
function. 

4. Since a is bounded below by —1 and decreasing the limit limr^oo o(r) = a exists and 
a G [—1, 0). For large r, and hence small /, Eq. ( p. 101) holds and implies that / is asymptotic 



to Cr where C is some positive constant. However, Eq. (|3.11| ) then tells us that 

for large r, a' ^ -{C^ /2\N\)r'^^~^\^\^''+'^^\ Thus a can only converge if 2|A^|(a + 1) - 1 > 1, 
or equivalently 

«>p-l. (3.12) 

This is impossible to satisfy if = 1, but for |A^| > 1 there is whole interval (— 1 + 1/|A^|, 0) 
of acceptable asymptotic values for a, and hence there is a one-parameter family of solutions 
of EVP, which was to be shown. 

When |A^| > 1 the variable parametrising the solutions of EVP can be taken to be a or, 
more physically, the total magnetic flux $ which is related to a via 

$ = J d'^xFu = N J d6a = 27iNa. (3.13) 

This formula shows in particular that the magnetic flux is not quantised. Mathematically 
this means that the gauge fleld {Ai, A2) of solutions of EVP does not extend to the compact- 
iflcation U {00}; for if it did, $/(27r) would be the flrst Chern number of a U{1) bundle 
over a compact manifold, which is necessarily an integer. The magnetic flux is nonetheless 
an interesting quantity to consider because it is conserved if one rules out inflnite energy 
conflgurations. This follows from Faraday's law of induction 



, E-dl, (3.14) 

dt Jc ^ ^ 

where the contour C is the circle at inflnity. The integral on the right hand side is only 
non-zero if the electric fleld E falls off for large r no faster than 1/r, which is precisely the 
condition for the electric fleld to have inflnite energy. 

We have numerically solved EVP with N = 2 for various values of $ in the allowed range 
(— 27r,0). Since the qualitative features of the functions / and a are clear from the proof 
of the proposition we only show plots of the energy density and the magnetic fleld of the 
solutions. The energy density is maximal on a ring whose radius is another useful measure 
of the soliton size. As the magnitude of the magnetic flux increases this radius increases and 
reaches a flnite limit for |$| 27i. In Figs. 1 and 2 the energy density and the magnetic 
field for a soliton whose (modulus of the) magnetic flux is close to that limit are drawn with 
a solid line. In the limit I'I'I ^ the soliton's size becomes arbitrarily small. 
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4 Discussion and outlook 



We have computed rotationally symmetric solutions of the Bogomol'nyi equations in a gauged 
version of the 0(3) sigma model of degree N with |A^| 7^ 1. There is no finite energy solution 
of degree 1, but there are probably many more solutions of degree > 1 than considered 
here. As in other field theories of Bogomol'nyi type this can presumably be shown using an 
index theorem and a vanishing theorem for an appropriate Dirac operator. Typically, there is 
a whole manifold of degree solutions of the Bogomol'nyi equation (called a moduli space) , 
and the dimension of this manifold is a linear function of A^. 

Solitons of Bogomol'nyi type which display all these properties and which are rather 
similar to the solitons discussed here are the Q-lumps mentioned in the introduction. Q- 
lumps of degree 1 necessarily have infinite energy, but there exists a (4A^ — 2)-dimensional 
family of Q-lumps of degree A" > 1. These include configurations which are made up of N 
well-separated single Q-lumps. Similarly there should be solutions of degree N > 1 in the 
present model whose energy density is peaked at A^ points in the plane. These could then 
be interpreted as superpositions of A^ solitons of degree 1; in this sense, solitons of degree 1 
can exist as part of a multisoliton configuration. 

Like Q-lumps the solitons discussed here can have an arbitrary size, with the role of the 
size parameter being played by the magnetic flux. However, whereas the energy of a Q-lump 
varies with the Q-lump's size, the energy of the solitons discussed here is degenerate with 
respect to changes in the magnetic flux. Thus the scaling degeneracy of the pure 0(3) sigma 
model persists in the present model in a mutated form. There is an important difference, 
however. In the 0(3) sigma model, solitons exhibit a "roUing instability", in the sense that 
under a small perturbation they either shrink to a thin spike or expand without limit (there 
is a subtlety here: in the so-called moduli space approximation such scale changes require 
inflnite energy for single solitons and most multisolitons; numerical simulations, however, 
suggest that they do occur in the full fleld theory). In the present model, by contrast, there 
is a completely general argument - Faraday's law of induction - according to which the total 
flux of a conflguration can only change at the cost of inflnite energy. 

The constancy of the total flux does not prevent individual solitons in a multisoliton 
conflguration from shrinking into thin spikes. It would be interesting to see if this happens 
as a consequence of soliton interactions. If the above conjecture about moduli spaces for 
solitons of degree N > 1 is correct this question could be investigated using the moduli space 
approximation to soliton dynamics. 

Finally we return to the question of stabilising solitons in the 0(4) sigma model in 
(3+1) dimensions with a Maxwell term. This would lead to the "Skyr me- Maxwell model 
without a Skyrme term" mentioned in the introduction. Of course it is well-known that 
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topologically non-trivial gauge fields can stabilise solitons in 2 or 3 spatial dimensions, like 
in the case of abelian Higgs vortices or t'Hooft-Polyakov magnetic monopoles. However, the 
the analysis of the gauged 0(3) sigma model here indicates that such a stabilisation is also 
possible if there is only a dynamical reason why the gauge field cannot vanish during time 
evolution. Thus one could try to stabilise a Skyrmion in three spatial dimensions via electric 
charge. The electric charge is conserved and, if it is non-zero, produces an electric field which 
prevents the Skyrmion from collapsing to zero size. It seems quite possible that there are 
such electrically charged Skyrmion solutions in (3-1-1 )-dimensional "Skyrme-Maxwell theory 
without at Skyrme term" . The basic idea is to balance the tendency of the scalar field to 
collapse by the electrostatic repulsion of like charges. 
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Figure Captions 



Fig. 1. The energy density as a function of r for solitons of degree N = 2 with $/(27r) = 
-0.96 (sohd line), $/(27r) = -0.43 (dashed hne) $/(27r) = -0.18 (dashed-dotted line). The 
function e plotted here is the integrand of ( p.7| ) divided by An. 

Fig. 2. The magnitude of the magnetic field F12 as a function of r for solitons of degree 
= 2 with <l>/(27r) = -0.96 (solid line), <l>/(27r) = -0.43 (dashed line) <l>/(27r) = -0.18 
(dashed-dotted line). 
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